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Abstract. The spectrum of temporal fluctuations of total magnetic energy for several 
dynamo models is different from white noise at frequencies smaller than the inverse 
of the turnover time of the underlying turbulent velocity field. Examples for this 
phenomenon are known from previous work and we add in this paper simulations of 
the G.O. Roberts dynamo and of convectively driven dynamos in rotating spherical 
shells. The appearance of colored noise in the magnetic energy is explained by simple 
phenomenological models. The Kolmogorov theory of turbulence is used to predict the 
spectrum of kinetic and magnetic energy fluctuations in the inertial range. 
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1. Introduction 

Several experiments have been carried out in the last decade in liquid sodium at 
high magnetic Reynolds numbers and in highly turbulent flows. Measurements of 
magnetic field fluctuations, either due to an externally imposed magnetic field or due to 
magnetic field generated through a dynamo effect by the sodium flow itself, generally 
reveal spectra with a l//-noise at low frequency [TJ [2j [3]. Numerical simulations of 
magnetohydrodynamic (MHD) flows showed the same phenomenon [H [5] although 
sometimes the precise form of the low frequency noise varies [6]. The precise value of the 
exponent in the algebraic decay in a spectrum is not so important, but a behavior other 
than white noise, i.e. colored noise, always invites closer inspection [TJ. For example, a 
1//- noise, frequently called flicker noise, cannot extend down to zero frequency if the 
power integral is to stay finite. The low frequency cutoff to the l//-spectrum should 
carry some informations about the dynamics of an MHD dynamo. 

Measurements of velocity in a fixed point of a turbulent flow frequently find white 
noise at low frequency (see for example the data compilation in fig. 6.14 of |8J), but there 
are also examples of colored noise [9]. The simultaneous presence of colored noise in 
both magnetic and velocity spectra looks like a plausible combination: If flow velocities 
increase in a dynamo, the magnetic field is amplified more rapidly and the amplitude of 
the magnetic field increases. If the velocity has no white noise, so should the magnetic 
field have a nontrivial spectrum. In the Karlsruhe experiment on the other hand [1] 
the l//-fluctuations in magnetic signals had no correspondence in pump pressures or 
volumetric flow rates. This prompts us to investigate in more detail the mechanisms 
leading to colored noise in magnetic spectra. 

This paper will deal with the fluctuations of total magnetic energy, which is a 
global measure of the field amplitude. A theory for the fluctuations of a global quantity 
is of interest even though laboratory experiments usually report local measurements. 
Numerical simulations always compute total magnetic and kinetic energies. Their 
statistics are less noisy than those of field amplitudes at a given point because of 
the spatial averaging inherent to the computation of total energies. Spectra of the 
geomagnetic field are frequently plotted as variations of another global quantity, the 
dipole moment pU] , because the geomagnetic field is dominated by the dipole component 
and low frequency variations reflect variations of the dipole moment. Secular variations 
of the Earth's magnetic field again exhibit colored noise [HJ [12] . In astrophysics, we do 
have local measurements of solar magnetic fields, but we only know a global amplitude 
for most stellar magnetic fields. 

The goal of this paper is to explore conditions leading to colored noise in the total 
magnetic energy and to show that this type of noise can appear even if the kinetic 
energy has white noise. Three simplified models of MHD flows are presented in section 
[2j The first two model the slow fluctuations per se, the last one by contrast computes the 
spectrum of inertial range fluctuations in homogeneous and isotropic turbulence. Section 
[3] presents numerical simulations of the G.O. Roberts dynamo [13J and compares the 



Global fluctuations in magnetohydrodynamic dynamos 



3 



results with section [2J whereas section H] does the same for simulations of convection 
driven dynamos in rotating spherical shells. 

2. Phenomenological models 

The basic problem in this paper is to find the temporal spectrum of a quantity derived 
from a magnetic field B(r,t), such as the magnetic energy. The induction equation 
governs the evolution of B in a liquid conductor with magnetic diffusivity A whose 
movement is given by the velocity field v(r,t): 



Exact solutions of this equation are difficult to obtain, so that we resort to 
phenomenological models. 

2.1. A single magnetic mode 

Mean field magnetohydrodynamics have proven a most fruitful simplification of the 
induction equation [TJ]. In this approach, the effect of small scale fluctuations on the 
large scales are not computed exactly but are modeled, in the simplest case as an 
a— effect. The number of magnetic degrees of freedom which need to be retained is 
thus reduced and in an extreme simplification, only one mode remains. If we call B 
the amplitude of that mode, a(t) and (3 the coefficients describing the a— effect and 
its quenching, respectively, and \i a coefficient related to the magnetic dissipation, the 
simplest model reproducing the main features of the induction equation is: 



a is allowed to be time dependent in order to reflect a time dependent velocity field. 
This time dependence will be essentially random for a turbulent velocity field. The 
reduction of the a— effect by the term (3B 3 models the retroaction of the magnetic field 
on the velocity field via the Lorentz force (which is quadratic in the magnetic field) in 
the Navier-Stokes equation. We now consider a(t) = a(t) — \x to be a random process 
with mean square (a 2 ) and remove the dimensions from eq. (T5]) by expressing time 
in multiples of (a 2 )" 1 and the magnetic field amplitude in multiples of ((a 2 )//?) 1 / 2 . 
The adimensional quantities t', a' and B' are given by t' = t(a 2 ), a' = a/ (a 2 ) and 
B' = B^J (3/ (a 2 ). In the remainder of this section, all quantities are understood to 
be nondimensional and the primes are omitted for convenience. The nondimensional 
variables then obey the equation: 



in which a(t) is a random variable with (a 2 ) = 1. As long as B is small, the solution to 
this equation is 




(1) 



d t B = a(t)B - (3B 3 - fxB. 



(2) 



d t B = a(t)B - B 3 



(3) 




(4) 
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For times small enough so that the exponent can be considered small, one has 



Taking the Fourier transform of this equation, it follows that the spectrum of B is, 
apart from frequency independent prefactors, the same as the spectrum of a divided by 
the square of the angular frequency, uj 2 . For example, if the spectrum of a is a white 
noise, the spectrum of B behaves as u~ 2 . For large times t, eqs. (jl]) and (E} become 
a poor approximation, which means that the u~~ 2 will not be observable below some 
cutoff-frequency. If the mean of a, (a), is different from zero, B will be large enough for 
the nonlinear term in eq. (|3]) to become dominant after a time on the order of (a) -1 . 
In that regime, and concentrating on slow fluctuations, eq. <^ reduces to 5 2 = a. 
Considering again the example of a(t) with a white noise, one finds a spectrum of B 
which is a white noise, too. The transition in the spectrum of B from u° to u~ 2 occurs 
at a frequency which increases with increasing (a), because eq. fl5]) fails at earlier times 
t. 

In order to test these ideas, we solved eq. ()3]) numerically. The random a(t) 
was generated by sending the output of a gaussian deviate random number generator 
through a Butterworth filter [151 [IB]- The filter was adjusted such that its output had 
a spectrum as a function of frequency / = uf (2tt) in 1/(1 + (/ / /i) 4 ) with fi = 50. The 
spectrum of B is shown in fig. [T]for different (a). As expected, this spectrum decays as 

for u> > u>i = 2nfi, and as u~ 2 for cu c < u < u\. Below the cut-off u c , the spectrum 
is independent of uj, and u c oc (a). 

Even though the spectra in fig. CD differ only in the value of u c , the qualitative 
appearence of the underlying time series is quite variable: The time series consists of 
intermittent bursts for small (a) and of random fluctuations around a well defined mean 
for large (a). A study of the probability distribution function of the solutions of eq. (J2j) 
together with a few graphs of representative time series is presented in [TF] . 

2.2. Several magnetic modes 

We will now investigate under which conditions the single mode model is applicable to 
more general systems and extend the discussion of the previous section to include several 
modes. It will be shown that the predictions of the single mode model are recovered in 
the limit of small fluctuations. 

The precise form of the dynamical system used as a model does not matter much 
for the following analysis, but a specific system has to be chosen for the numerical 
examples. In order to stay as close as possible to the previous section, let us assume 
V • v = and rewrite the left hand side of eq. ([T]) as d t Bi + J2j v jdjBi — J2j BjdjVi. We 
then proceed through the same steps as before and replace the combination of velocity 
and derivation by a random variable in which we absorb the dissipative term and the 
right hand side, remove dimensions, and model saturation through a cubic term. This 




(5) 
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leads to the following system: 

d t B i +(a 1 {t)+a 2 (t)+a 3 {t))B i -a i (t)(B 1 +B 2 +B 3 ) = -Bf , t = 1,2,3, (6) 

in which the a,(t) are random variables. This system bears only a metaphorical relation 
with the original induction equation and will be used to exemplify two limiting cases: 

If the fluctuations of the ctj(t) are small compared with the mean of the the 
solution of (EJ) will be close to the solution of the time independent system in which each 
o>i(t) in ([6]) is replaced by its mean (a,). Let us assume = (a 2 ) — (0*3) < 0. An 
eigenvalue analysis of the left hand side of (EJ) then reveals one neutral mode and two 
modes with equal and positive growth rate. In the presence of small fluctuations, the 
neutral mode will not contribute significantly to the dynamics. If the initial conditions 
and the nonlinear term select an arbitrary direction in the space spanned by the two 
degenerate growing modes, we expect ([6]) to behave the same as the single mode model. 
This is a fortiori true for a dynamical system with a single non-degenerate growing 
mode. 

If on the other hand the fluctuations of the are large compared with their 

means, the dynamics is not dominated by a single mode anymore and the analysis of 
the previous section breaks down. The spectrum of the fluctuations of \Jj2 Bf can now 
be different and must be found from numerical computation. 

Fig. [2] shows some examples of solutions of eq. ([6]) in which the spectrum of the 
fluctuations of the ai(t) is in 1/ / [18] and (a?) = 1. Let us first consider the case in 
which the fluctuations of the a, are small compared with their mean. The spectrum of 
\Jj2 Bf shown in fig. [2] should then behave as predicted by the single mode model: At 

the smallest frequencies, the spectrum of \Jj2 Bf must decay the same as the spectrum 
of the ctj, i.e. as 1/f in the present example. Above a frequency on the order of 
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there must be a factor f 2 between the power laws followed by the spectra of Jj2 Bf and 

the CKj, which implies a spectrum in 1/ f 3 for \Jj2 Bf in the example considered here. All 
these predictions fit well the spectrum shown in fig. [2] for (a 2 ) = 1 and (c^) = —5. In 
order to further support the applicability of the single mode model, we also computed 
the angle between the instantaneous vector B(t) = (Bi(t) , B 2 (t) , B 3 (t)) and its mean 
(B). The cosine of that angle, (B) • B{t) / ^\{B)\ 2 \B(t)\ 2 in the statistically stationary 
state stays larger than 0.99 for (a^ = —5 but is scattered over a large interval for 
(a^) = —0.01. In the latter case the fluctuations of the are large compared with their 
mean and different exponents unrelated to the single mode model become possible. A 
decay in 1/f 2 appears in fig. [2] which will become important again in connection with 
the convectively driven dynamos in spherical shells discussed below. 




f 



Figure 2. Spectral power density of Bf obtained from the solution of cq. ([5]), as 
a function of frequency / = lj/(2tt), for (af) = 1 and (ai) = —0.01 (red), -0.3 (green) 
and -5 (blue). The straight lines indicate the power laws to , co~ 2 and u>~ 3 . 



2.3. Homogeneous and isotropic turbulence 

The main focus of this paper is on slow fluctuations, but it is also worthwhile to have 
a look at the fast fluctuations in order to see where the differences are. We also would 
like to check whether it was reasonable in section [2TT1 to assume a spectrum for a which 
is flat at low frequencies and steeper at large frequencies. We will specifically look at 
fluctuations of total kinetic and magnetic energy. The spectra of these fluctuations 
must be different from flow to flow, but we can expect a unique behavior in flows to 
which the Kolmogorov phenomenology (K41 in short after [19]) of homogeneous and 
isotropic turbulence applies. We will employ this phenomenology in a cartesian domain 
with periodic boundaries, so that it is useful to expand all fields in Fourier series. Much 
effort has been spent in the past to theoretically deduce the wavenumber dependence 
of the Fourier coefficients, the best known theory being of course K41. But the Fourier 
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coefficients also depend on time, or after an additional transformation, on frequency. 
Theories such as K41 are concerned with the temporal mean of those coefficients. The 
spectrum of temporal fluctuations of those Fourier coefficients has received much less 
attention. 

This section will present an extension of the Kolmogorov phenomenology which 
predicts the spectrum of fluctuations of the total kinetic energy of a flow. There is no 
magnetic field involved, but the result is directly applicable to magnetic field spectra 
under certain assumptions. 

The mean kinetic energy (Skin) of a turbulent flow is decomposed into contributions 
made by wavevectors of modulus k in the form (£kin) = Jo°(E k )dk. According to K41, 
(E k ) oc A; -5 / 3 . The number of modes in a thin shell of radius k in spectral space is 
proportional to k 2 , so that the contribution of a single mode with wavevector k to 
the total mean kinetic energy is proportional to k^ 11 ^ 3 . We now consider the root 
mean square of the fluctuations of the total kinetic energy: a tot = ((E^ in — (E^) 2 )) 1 ^ 2 . 
The angular brackets denote average over time. The total root mean square is again 
decomposed into contributions of different wavevector shells: a to t = Jo°° °~kdk. In order 
to compute a k , we assume that the fluctuations of the k 2 modes in the wavevector shell 
of radius k are all uncorrected so that the mean squares of the fluctuations of each 
single mode simply add to give a\. 

We next invoke the concept of self-similarity, which is a central tenant of the 
K41 theory: All structures (or eddies or modes) in the inertial range are statistically 
indistinguishable from each other after a rescaling of length and time. This implies that 
the histogram of the fluctuations in every mode has the same shape if the fluctuations are 
given as multiples of the mean amplitude. We then have to conclude that the root mean 
square of the fluctuations scales the same as the mean amplitude, i.e. is proportional to 
k~ 11/s . It follows that a k oc VWk' 11 / 3 = k~ 8/3 . 

Note that it is not possible to deduce the scaling of a k from the usual dimensional 
arguments of K41. The ratio a k j (E k ) tends to zero if the number of modes in a 
wavevector shell tends to infinity because the fluctuations of the modes in that shell 
average out to zero. The rms of the total fluctuation depends therefore on the number 
of degrees of freedom, which in turn depends on the ratio of integral to dissipative length 
scale. These two length scales are assumed to be irrelevant in the K41 theory, however. 

The energy of the modes in the wavevector shell of radius k is now written in the 

form 

E k (t) = (E k ) + a k h k (t) (7) 

in which h k (t) has to obey (h k ) = and (h 2 k ) = 1. The Fourier transform of the 
fluctuations of the total energy is 

J (E k (t) - (E k ))e-^dt oc J a k h k {uj)dk (8) 

with h k (uo) oc / h k {t)e~ %UJt dt. We now use again the hypothesis of self-similarity and note 
that the typical time scale for a mode with wavenumber k is its turn-over time which is 
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proportional to k 2 / 3 [2TJ] IS]. Restricting ourselves to frequencies u> and wavenumbers fc 
in the inertial range, we expect 

h k (u) oc k~ 1/3 g(ouk- 2/3 ) (9) 

with an unknown but universal function g. The form of the argument reflects that 
all hk{ui) should have the same form once u is expressed in multiples of the turn-over 
frequency, and the amplitude factor k" 1 ^ 3 follows from the requirement that (h 2 ,) = 1, 
which implies that / |/ifc(c<;)| 2 <ic<; is a constant independent of k. Substituting z = uk~ 2 ^ 3 
into J \h k (u)\ 2 du = J k~ 2 / 3 \g(uk~ 2 / 3 \ 2 )duj shows that this is indeed fulfilled. 

We can now proceed to finally compute the spectrum of the energy fluctuations, 
using the same substitution z = cuk~ 2 / 3 , to find: 

J {E k (t)-(E k ))e- iult dt oc j k- % l 3 k- l ' 3 g(uk' 2 l 3 )dk oc u~ 3 J z 2 g{z)dz.(10) 

The spectrum of the energy fluctuations, i.e. the square of the Fourier spectrum above, 
must thus decay as w~ 6 for u in the inertial range. This prediction will be verified below. 

Spectra in MHD turbulence can have a variety of shapes, depending on the presence 
of Alfven waves, an applied magnetic field, the magnetic Prandtl number etc. [2T| 122] 
The above calculation can be directly reproduced for the magnetic energy as long as 
the magnetic energy density follows a Kolmogorov spectrum in k~ 5 ^ 3 . This happens in 
homogenous and isotropic flows at high magnetic Reynolds numbers without applied 
external field [23]. In that case, magnetic field fluctuations should decay as w -6 , too. 



3. The G.O. Roberts dynamo 



Dynamos based on periodic flows first investigated by G.O. Roberts [13J have been 
useful for a number of basic studies of the dynamo effect [211 [25j [261 123 |28j and have 
inspired the Karlsruhe dynamo experiment [29]. In order to verify the ideas developed 
in the previous section, we wish to numerically simulate the Navier-Stokes and induction 
equations in the following non-dimensional form: 

Pm -°" (11) 

(12) 

(13) 



d t B + V x (B x u 
V-B = 
d t u + u ■ Vu 



Re 



-V 2 B 



-Vp + t-V 2 u + F + (V x B) x B 
Re 



V ■ u = 



(14) 



Re and Pm are two control parameters standing for Reynolds and magnetic Prandtl 
numbers, respectively. The forcing F will take the form 
_ 8vr 2 



Re 



■Uq 



One solution to the above equations is then B 

^ \pl sin(27rx) cos(27n/) \ 

— y/2 cos(27rx) sin(27n/) 
, 2 sin(27ra;) sin(27r?/) j 



0, u = Uq with 



Uq 



(15) 



(16) 
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This solution is unstable for large enough Re, either because the flow is 
hydro dynamically unstable or because the dynamo effect leads to B / 0. The velocity 
field ( fl6l) is a cartesian arrangement of helical eddies with their axes along z. Periodic 
boundary conditions will be applied in all cartesian directions x, y, z with periodicity 
lengths l x = 1, l y = 1, and three different choices of l z , namely l z = 1, 1.5 and 2. 

If we call in dimensional units L x the periodicity length in x— direction, the laminar 
velocity v = VoU , and the kinematic viscosity and magnetic diffusivity v and A, 
respectively, the nondimensional parameters are given by Pm = vj\ and Re = VqL x /v. 

The above problem is conveniently discretized with a spectral method. However, 
we need very long time series if we are interested in small frequencies and the availability 
of computer hardware for long runs becomes an important consideration in the choice 
of the numerical method. The most readily available platform happened to be a set of 
machines equipped with CUDA (compute unified device architecture) capable graphic 
processing units. These architectures are not well adapted to spectral methods, nor to 
Poisson solvers. That is why a fully explicit finite difference scheme was implemented 
to solve the following set of equation: 

r m 

d t B + V x (B x u) = —V 2 B (17) 

Re 

V-B = (18) 

d t u + u -Vu = -c 2 Vp + ^-V 2 n + F + (V x B) x B (19) 

Re 

d t p + V • u = (20) 

These equations describe a compressible fluid with the equation of state p = c 2 p and a 
linearized continuity equation. Solutions of the incompressible Navier-Stokes equation 
will be recovered in the limit of the sound speed c tending to infinity. The sound speed 
was chosen large enough so that the peak Mach number never exceeded 0.12 in the 
computations presented below. The typical Mach number was more like 0.05. The 
finite difference code was validated against a spectral code [30] by comparing the results 
of short runs, but all the heavy computation was done with the finite difference code. 

We will first look at non-magnetic flows. Figure |3] shows spectra of kinetic energy 
density E^ n for Z 2 = 2, 1.5 and 1. E\^ n is defined as 

1 r , 



£kin = J u 2 dV, (21) 

where V is the computational volume. In all cases, the spectra of E^ in follow an 
a; -6 over at least a decade in u>. This power law identifies the inertial range of a 
Kolmogorov cascade according to section 12.31 (It was not checked directly whether the 
spatial spectrum of the velocity fields follows the K41 law because the simulations were 
not done with a spectral code and the evaluation of the spatial spectra would have been 
cumbersome). An even steeper decay follows at higher frequencies before the spectrum 
dives under the noise level introduced by roundoff error. 

On the low frequency side, the inertial range ends in a conspicuous hump. The 
frequency of the local maximum in the spectrum, / max , scales with what one may identify 
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Figure 3. Spectral power density of the kinetic energy as a function of frequency 
/ = uj/(2tt) for flows without magnetic field. The panels show spectra for l z = 2, 1.5 
and 1 and each panel contains spectra for Re = 10 3 (red), 3 x 10 3 (green), 10 4 (blue) 
and 3 x 10 4 (pink). The straight lines indicate the power laws w" 1 and w~ 6 . 
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Table 1. Z z , Re, Re2 and Re2/ m ax for the simulations in fig. [3] 

as the energy injection scale: Re is a control parameter in eq. (fT4l) . but the Reynolds 
number determined a posteriori is Re 2 = V-^kin- The values of Re 2 given in tabled] vary 
from 340 to 2100 indicating that all flows are turbulent. The product Re 2 / max varies 
from 27 to 41. Despite this variation, we identify / max with the injection scale because 
the definition of Re 2 does not take into account that the flow is anisotropic so that we 
do not expect Re 2 / max to be strictly constant. 

The spectra at frequencies smaller than the injection frequency are different for 
l z = 1 and l z = 2 or 1.5: For l z — 1, the spectrum is purely white noise, whereas a 
section of l//-noise is visible for l z = 2 and 1.5. The section of l//-noise shrinks with 
increasing Re. 

We next turn to dynamos. For Pm slightly above the critical value given in table 
[U the spectra of the magnetic energy density E B , defined as E B — ^ J B 2 dV, show a 
segment of spectrum in u~ 3 for l z = 2 and 1.5 and in u~ 2 for l z — 1. A white spectrum 
is always found at the smallest frequencies (figure H]). This matches the prediction of 
section [27T1 in as far as the velocity spectrum without magnetic field contains an interval 
with a spectrum close to u^ 1 for l z — 2 and 1.5 and nothing but white noise below the 
injection frequency for l z — 1. According to section 12. 1[ there should be a factor u 2 
between magnetic and velocity spectra at low frequencies, which is compatible with the 
results in figures [3] and HI 

The transition from white noise to either u~ 2 or u~ 3 occurs at a higher frequency 
for larger Pm in figure HI This is again in agreement with section 12.11 because a larger 
Pm at constant Re corresponds to a flow farther above onset, which corresponds to a 
larger (a) in the model of section 12.11 

Experiments use local probes to characterize the magnetic field. In the present 
simulations, the magnitude \B\ of the field at position x — y — 0, z — l z /2 was recorded 
as a function of time. This location corresponds to a position at which the magnetic 



Global fluctuations in magnetohydrodynamic dynamos 



12 




f 



Figure 4. Spectral power density of the magnetic energy Eb as a function of frequency 
/ = u)/(2w). The panels show spectra for l z = 1.5 (top) and 1 (bottom) and Re = 10 4 . 
The different lines in the top are for Pm = 0.03 (red), 0.05 (green) and 0.1 (blue), and 
in the bottom panel for Pm = 0.27 (red), 0.5 (green) and 1 (blue). The straight lines 
indicate the power laws oj~ 3 and co~ 6 (top) and oj~ 2 and w -6 (bottom). 

field has been measured in the Karlsruhe experiment pQ. The spectra of the magnetic 
field amplitude in one point show a power law at low frequencies over a larger interval 
than the spectra of Eb, but with an exponent which bears no obvious relation with the 
exponent of the energy spectra (figure [5]). There does not seem to be any theoretical 
tool for predicting local fluctuations, so that we simply note the spectra in figure [5] as 
an empirical fact. However, they also approach a white noise spectrum as Pm increases, 
just as the magnetic energy spectra do in fig. HI 

4. Dynamos in spherical shells 

Simulations of convectively driven dynamos in spherical shells reported in [6] yielded 
spectra with power laws which are not compatible with the single mode model. 
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Figure 5. Spectral power density of the local magnetic field amplitude \B\ at the 
position x = y = 0, z = l z /2 as a function of frequency / = u>/(2tt). The panels show 
spectra for l z = 1.5 (top) and 1 (bottom) with Re = 10 4 and the same Pm as in fig. 
HI The straight lines indicate the power laws (top) and w" 1 / 2 (bottom). 

According to the phenomenology developed in section 12. 2\ this simply means that the 
fluctuations in these dynamos are large. It is expected that at sufficiently low Rayleigh 
numbers and in sufficiently quiet convection, one finds again spectra in agreement with 
the single mode model. In order to test this idea, we also performed simulations of 
dynamos in spherical shells. Because these computations are more expensive, the spectra 
are more noisy and not as extended in frequency as for the G.O. Roberts dynamos, but 
they demonstrate the main effect. 

We consider the same physical model as in [6] of a rotating spherical shell with its 
gap filled with fluid driven by convection. The variable quantifying buoyancy will be 
called temperature here, but the employed boundary conditions better fit compositional 
convection: Fixed temperature at the inner boundary and zero heat flux through the 
outer boundary. Both boundaries are assumed no slip and electrically insulating. The 
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following dimensionless equations 

E 2 Ra r 

d t u + (u ■ V)u + 2zxu = V<f) T— + (V x B) x B + EV 2 u(22) 

Pr r 

d t T + (u ■ V)T = ^V 2 T - Ee (23) 
Pr 

d t B + V x (B x u) = — V 2 B (24) 

Pm 

V-w = , V ■ B = (25) 

are solved in a spherical shell for the velocity, magnetic and temperature fields u, B 
and T. The control parameters are the Ekman number E, the Rayleigh number Ra, the 
Prandtl number Pr, and the magnetic Prandtl number Pm. The ratio of inner and outer 
radii is fixed at 0.35 with the dimensionless outer radius being tq = 1/0.65, and e, the 
variable modeling the buoyancy source [6] is fixed at e = 1. The equations have been 
solved using the spectral method described in [3]] with a resolution of 33 Chebychev 
polynomials in radius and spherical harmonics of degree up to 64. 

Fig. [6] shows two calculations in which all parameters are held constant except for 
the Rayleigh number. In going from the low to the high Rayleigh number, the spectrum 
of kinetic energy changes from white noise to approximately The spectrum of 

magnetic energy shows a decay in f~ 2 in both cases. This means that at low Ra, the 
exponents in the power laws for the kinetic and magnetic energies differ by 2, just as 
they should according to the single mode model. At higher Ra, i.e. in the flow with 
larger fluctuations, the exponents are different as expected from section 12.21 and differ 
by 1. Incidentally, the combination of f^ 1 and f~ 2 for kinetic and magnetic energies is 
identical with a combination of exponents found in the low dimensional dynamic system 
of section 12.21 



5. Conclusion 



Colored noise has been noted in fluctuations of magnetic field strength in both 
experiments and simulations. According to common experience with turbulent velocity 
fields, the fluctuations of mechanical quantities have on the contrary a white noise at low 
frequencies, a low frequency being a frequency smaller than the inverse of the integral 
time scale. It has been shown in the present paper that even the simplest phenomenology 
based on a single mode model predicts a spectrum in 1/ f 2 for the total magnetic energy 
if the velocity field is characterized by white noise. The single mode model is justified 
as long as the fluctuations in the velocity field are small enough in amplitude. These 
simple models are corroborated by numerical solutions of the dynamo equations for 2D 
periodic dynamos and for convectively driven dynamos in spherical shells. 

Colored noise must have a low frequency cut-off for the power integral to converge. 
This low frequency is independent of diffusive time scales according to the models studied 
here. Instead, it depends on how strongly a dynamo is driven. The more supercritical 
a dynamo is, the larger is the cut-off frequency. 
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Figure 6. Spectral power density of the kinetic energy E^m (top) and the magnetic 
energy Eb (bottom) as a function of frequency / = u>/(2tt) for Pm = 13, Ek = 
6.5 x 10~ 3 , Pr = 1 and Ra = 10 5 (red) or Ra = 10 6 (green). The straight lines indicate 
the power laws /° and f" 1 in the top panel and f~ 2 in the bottom panel. 

It is also possible to deduce the spectrum of turbulent magnetic fluctuations 
from the Kolmogorov phenomenology. Within the inertial range, fluctuations of the 
total kinetic and magnetic energy should decay as f~ 6 , which is verified by numerical 
simulation. 

It is left to future work to deal with fluctuations of the magnetic field at a given 
point in space. If the magnetic field is dominated by a large scale component, local 
fluctuations are of course similar to global fluctuations. At large magnetic Reynolds 
numbers, small scale fluctuations exceed in magnitude the fluctuations of the large scale 
field and the spectrum of the fluctuations of local magnetic fields can be different from 
the spectrum of total magnetic energy. The spectra computed for the 2D periodic 
dynamos at high magnetic Reynolds numbers reveal power laws at low frequencies in 
the range to /~ ' 5 . 
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